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Abstract 

A dynamical decoupling scheme for the deterrence of errors in the non-Markovian 
(usually corresponding to low temperature, short time, and strong coupling) regimes 
suitable for qubits constructed out of a multilevel structure is studied. We use the 
effective spin-boson model (ESBM) introduced recently [K. Shiokawa and B. L. Hu, 
Phys. Rev. A70, 062106 (2004)] as a low temperature limit of the quantum Brownian 
oscillator model, where one can obtain exact solutions for a general environment with 
colored noises. In our decoupling scheme a train of pairs of strong pulses are used to 
evolve the interaction Hamiltonian instantaneously. Using this scheme we show that 
the dynamical decoupling method can suppress 1// noise with slower and hence more 
accessible pulses than previously studied, but it still fails to decouple super-Ohmic 
types of environments. 
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1 Introduction 



In quantum information science, quantum coherence and entanglement are used as essential 
resources for efficient information processing. Since the interaction of a system with its 
environment can destruct the coherence and destroy the entanglement, it is regarded as 
the most serious obstacle in the realistic implementation of quantum information processing 

Fast time low temperature strong field challenges In our previous paper |4;, we 
introduced an Effective Spin-Boson Model (ESBM), obtained as the low temperature limit 
of the quantum Brownian oscillator model (QBM) 0101 when the system behaves effectively 
like a two-level system. In the absence of external fields, the solutions of ESBM match 
with the exact solutions for the spin-boson model[71 |H]. At a finite temperature and in 
the presence of an external field, the multilevel structure inherent in ESBM provides a 
quantitative measure of the leakage from the two-level system by transitions to other levels. 
Although leakage is unavoidable in many realistic qubit systems, most theoretical models 
ignore it. It is either interpreted as the difference between the theoretical prediction and 
the observation or estimated phenomenologically by crude approximations such as invoking 
the Fermi Golden rule which makes the Markovian assumption which may contradict the 
existent conditions. 

Our model based on exact solutions for QBM enables us to estimate leakage analytically 
and provides an essential tool for probing the multilevel coherent dynamics without invoking 
the conventional approximations such as Born, Born-Markov, rotating-wave approximations 
(RWA). It also provides a systematic method for studying the quantum coherence of the 
system behavior in the presence of arbitrary external pulses which is used in error correction 
schemes. The exact decoherence and Rabi-oscillations in Ohmic and super-Ohmic environ- 
ment were studied in j3]. We emphasize that we don't invoke the Markov approximation 
which is likely to be unreliable in the presently studied qubit models operating at low tem- 
peratures. In order to study the short time dynamics under strong and fast pulses, RWA, 
obtained by ignoring the counter rotating terms averaged out at long times, is also not suit- 
able. Our model does not depend on these approximations and can thus serve the special 
yet important purpose of analyzing the short time dynamics under strong external fields. 

Error Correction versus Deterrence One of the proposed methods to remove errors 
from quantum computation architecture is the use of error correcting codes, which enable one 
to detect and correct errors after they have occurred|H]. By encoding qubits in a redundant 
way, we can detect errors from majority voting. Although this method works irrespective of 
the source of errors, actual coding requires many physical qubits including ancillas and its 
implementation is not practical at the present stage of solid state qubits. When the system- 
bath interaction occurs collectively, we can encode qubits in the decoherence free subspace 
in order to avoid decoherence before it occurs|H|. 

In this paper, we apply the ESBM to quantum error deterrence by dynamical decoupling. 
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in which the external pulses are used to overpower the influence of the environment on the 
system. Using the ESBM, we give a more systematic study of the decoupling method which 
hopefully provides deeper insights in these processes than previously done based on simpler 
models. 

In the popular scheme known as "bang-bang" (BB) control one introduces a sequence of 
square pulses [inilIIlliaiISllIlllI3linilIZllIHl- Most theoretical analysis of this method are 
directed at pure-dephasing in the spin-boson model: the system-bath interaction is through 
the Pauli matrix cx^, or the atom-field interaction under the rotating wave approximation [13] . 
Similar techniques are used for the study of suppression of collisional decoherence [20] and 
spontaneous decay (TH] or the modified coherent dynamics by applying sinusoidal waves pTj. 

Dynamical decoupling by spin-echo pulses has been commonly used to suppress deco- 
herence due to magnetic field fluctuations for ion-trap quantum teleportation Similar 
to the spin-echo a pair of vr pulses are used in the bang-bang control. In contrast, our de- 
coupling method uses a rapid sequence of a pair of pulses with opposite sign to rotate the 
oscillator coordinate in the same direction in the phase space by mimicking the instanta- 
neous Hamiltonian evolution. Since our ESBM model has more realistic features than those 
models previously used to study decoupling and the advantage of being exactly solvable, it 
provides an ideal testing ground for the effectiveness of dynamical decoupling in realistic 
settings beyond the regimes delimited by conventional approximations. 

It has been argued that for the decoupling to be successful one needs to use fast pulses 
with pulse duration much shorter than the inverse frequency cutoff of the bath. 



In |2S1, it was shown that when dynamical decoupling is applied to 1// noise, or more 
generally, 1//° (a > 0), the pulses can be somewhat slower although the condition 



still needs to be satisfied. The decoupling was shown to be more effective when the noise 
originates from a bath with 1// spectrum than in the Ohmic case due to infrared (IR) 
dominance of modes in the 1// spectrum. 

Noise with 1// type of spectrum are ubiquitous and appears in a wide variety of cir- 
cumstances, for instance, river flow, sunspot activity, neuronal spike trains, and human 
coordination^^. Various interpretations are given but so far no explanation is universal. 
For electronic systems, it is often attributed to the motion of defects or impurities or back- 
ground charge fluctuations, although different interpretations are possible. 

In j2Hl! an experiment was performed with spin-echo- type pulses to identify the dominant 
source of decoherence in superconducting charge qubits made of Cooper-pair boxes. The 
authors concluded that a dominant source of decoherence is due to 1// charge noise. In 
the authors argued that similar techniques can be used to efficiently suppress the 1// noise. 
They also discussed the limitation and general applicability of decoupling pulses. Similar 
results were found in the study of bang-bang control against telegraph noise |25|. 




At < 1/Ac/y, 




At < n/Auv, 



(2) 
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On a short time scale comparable to the bath time scale, the details of the system-bath 
interaction and internal bath dynamics are no longer negligible. Those details will emerge 
when non-Markovian features of the dynamics become important. For the QBM model, these 
details are captured by the bath spectral density /(a;). For 1// noise the spectral density 
is concentrated in the low frequency range. As a result, the decoupling result depends 
more sensitively on the lower than on the upper cutoff. For the pulse interval close to the 
threshold value, decoherence rate can increase (see Eq. (4)). Similar accentuation of 
the decay process can occur from infrequent measurements, which is known as the quantum 
anti-Zeno effect |26i|. 

In Sec. 2 we discuss the effect of a quantum Brownian harmonic oscillator (QBO) under 
the influence of a general environment and an external pulse field. We use the Fock space 
representation for the harmonic QBO to define an effective two level system. In Sec. 3 
we introduce dynamical decoupling (DD) in this effective spin-boson model (ESBM) and 
describe how to design pulses to minimize the influence of the environment and limit the 
leakage (to higher levels) . We describe how well our DD scheme work in error deterrence for 
different environments. In Sec. 4 we draw the conclusions. 



2 Modified quantum Brownian motion due to pulses 
2.1 The model 



Our model consists of a quantum Brownian particle interacting with a thermal bath in the 
presence of external pulses. We use the units in which = h = 1. The Hamiltonian for 
our model can be written as 



where the dynamics of the system S (with coordinate x and momentum p) is described by 
the Hamiltonian 



The Hamiltonian of the environment is assumed to be composed of harmonic oscillators with 
natural frequencies a;„ and masses m„. 



where {qi, ...,gjv,Pi, ■■■,Pn) are the coordinates and their conjugate momenta. The interac- 
tion between the system S and the environment E is assumed to be bilinear. 



H = Hs + He + Hj + Hp, 



(3) 




(4) 




(5) 



N 



Hi = Cnq, 



(6) 



n=l 
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where c„ is the couphng constant between the quantum Brownian oscillator (QBO) and the 
nth environment oscillator with coordinate g„. The coupling constants enter in the spectral 
density function I{uj) of the environment defined by, 



JM^vr5:-^5(a;-^„). (7) 

We assume the spectral density has the form 

/(cu) =2M7cu"e-"/^^^, (8) 

where u = 1 is Ohmic, z/ < 1 is sub-Ohmic, and u > 1 is supra-Ohmic. We will study 
Ohmic, super-Ohmic with u = 3, and sub-Ohmic with 1/f type noise. Since the last case 
is sensitive to the infrared (IR) regime, we introduce an infrared (IR) cutoff Ajr for our 
frequency integral. The law of large numbers tells us that a general environment consisting 
of a large number of fluctuators weakly coupled to the system can also be reduced to a 
bosonic environment. This description is supported numerically and has been applied 
successfully to explain the effect of 1// noise on superconducting charge qubitsj2H]- 
For a linear QBO, the potential is 

Vi^) = , (9) 

where fl is the natural frequency of the system oscillator. The pulse Hamiltonian has the 
form 

N 

Hp = HsY.Pn{t) (10) 

n=l 

where 

At <t < At + T (mod 2 At + 2r) 
P„(t) = \ -V -T <t<0 (mod 2At + 2r) 

otherwise 

Thus T is the duration of the pulse and At is the pulse interval between the end of one pulse 
and the beginning of another of opposite sign. 




2.2 The influence functional of a system acted on by pulses 

Here we will make connection with prior treatments of QBM based on the influence functionalj3] 
with a phase space representation for the Wigner function |30J. First we define the transition 
amplitude between the initial state \xo Qq) a.t t = and the final state \x q) at time t to be 

K{x,q;t I Xo,qo;0) = {x q\U{t,0)\xo go)- (H) 
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where U {t, 0) = e~*^* is a time evolution operator. In the presence of strong pulses, i.e. V — > 
oo while keeping Vt fixed, the pulse Hamiltonian commutes with the system-environment 
evolution Hamiltonian. Then the total amplitude can be sequenced in time as 

K{x,q;t I xo,qo;0) = {x q\U{t,t2N)UlU{t2N,t2N-i)UpU{t2N-i,t2N-2)-- 

UlU{t2,h)UpU{h,0)\xo qo) (12) 

= {x q\U{t,t2N)Up{t2N,hN-l)U{t2N-l,hN-2) — 

Up{t2,ti)U{ti,0)\xoqo) (13) 

where Up = e''^^^^ and Up{t2N,t2N~i) = e'"s^^U{t2N ,t2N-i)e-'"^^^ is the modified evo- 
lution kernel due to pulses. This form shows that the effect of a pair of pulses is to 
evolve the system variable x in the interaction Hamiltonian instantaneously to xcos{VtVT) + 
psin(QyT)/Q. In phase space, this corresponds to an instantaneous rotation by an angle 
VtVr. 

The Liouville equation for the density matrix is 

^^p{t) = [H,p{t)], (14) 

where [, ] is the commutator. In the coordinate representation, the density matrix becomes 

p{x, x\ q, q', t) = {x q\p{t) \x' q') (15) 

where we have used q = {q^} to denote the environment variables collectively. The time 
evolution of the density matrix is given by 

p{x, x', q, q, t) = j dxQdxQdqod%K{x, q; t \ Xo, qo; O)p{xo, x'q, qo, q'^, 0)K*{x', q; t \ x'q, q'^; 0).(16) 

For a factorized initial condition between the system and the environment: p{xo, Xq, qo, qQ, 0) = 
Ps{xo,Xq,0) <S> PB{qo,(lo,0), after tracing out the environment harmonic oscillator variables, 
we obtain an equation for the reduced density matrix: 



Pr{x, x', t) = J dqp{x, x', q, <l,t) = J dxodx'^Jrix, x'; t \ Xq, x'^; O)ps{xo, x'^, 0), (17) 
where its time evolution operator is given by 
Jr{x, x'\ t I Xo, Xq; 0) = y dqdqodqQK{x, q; t \ Xo, qo; 0)ps(go, ^o, 0)K*{x', q; t \ Xq, q'^; 0). (18) 



The total action of the system S[x, x'] enters as 

r{xx') 



Jr{x, x'; t I Xo, 0) = / PxPx'e^'^l^'^ 1 

(19) 
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and consists of two contributions: 

S[x,x'] = Ss[x,x']+Sif[x,x']. (20) 
The action for the system Ss is given by 

Ss = /* ds{MR(s)f(s) - Mn^R(s)r(s)}, (21) 
Jo 

where R = {x + x')/2, r = x — x'. The influence action Sif[x,x'] accounts for the effect of 
the environment on the system S and is given by 

SiF[R,r]—i [ ds [ ds'y{s)i'{s — s')y{s') 
Jo Jo 

-2 [ ds r ds'y(s)i^(s- s')Y(s'), (22) 
Jo Jo 

where y(t) = {x{t)+x'{t)) cos 9{t)/2+M{x{t)+x'{t)) sin e{t)/2n, y{t) = {x{t)-x'{t)) cos9{t)+ 
M{x{t) — x'{t)) sm9{t)Q. The time dependent angle 9{t) is defined as 

n(.. _ / At<t<2At (mod 2At) , . 

'^^^^ ~\QVt -At<t<0 (mod 2At) ^ "^^ 

and 

1 f°° 

fiU) = / du>I(uj)smujt (24) 

TT JO 

vit) — — I duliijj) cot\i — cosut, (25) 

TT JQ 2 

are the dissipation and noise kernels respectively for an environment initially at thermal 
equilibrium. We are particularly interested in the case in which each pair of pulses rotate 
the system by an angle tt in the phase space, i.e., QVt = {2n + l)7r with an integer n. In 
this special case, due to the periodic nature of the harmonic motion, difference between our 
method and the bang-bang control is reduced to the phase factor. The effect of these pulses 
is to modify the form of the influence kernels in the action to 

SjF[R,r]=i ds ds'r(s)i>(s,s')r(s') 
Jo Jo 

-2 t ds r ds'r{s)fl{s,s')R{s'), (26) 

Jo Jo 

where the modifled kernels jl and i> are deflned as 

il(s,s') = cos9(s)ii(s-s')cos9(s') (27) 
i>{s,s') = cos9(s)u(s- s') cos 9(s') (28) 
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^From Eqs. ()26p the Euler-Lagrange equations for R and r are 



MR(t) + MQ^R{t) + 2 /* ds^{t - s)R{s) = 0, (29) 

Jo 



t 



Mr{t) + Mn'r{t) - 2 / ds'j2{t - s)r{s) = 0. (30) 

Js 

The effect of the environment is included in the nonlocal kernels /i and u. On the other hand, 
these kernels contain the information of the past history of the environment as modified by 
the presence of the system variables. Thus these equations respect the self-consistency in 
the dynamics between the system and the environment. 

If we formally write the two independent solutions of the homogeneous parts of Eq. pHj) 
and Eq. (jHUjl to be Ui{s) and fj(s), i = 1,2, with boundary conditions Mi(0) = l,Ui(t) = 0, 
^2(0) = 0,U2{t) = 1 and fi(0) = l,Vi{t) = 0, ^2(0) = 0,V2(t) = 1, respectively, the solutions 
of these uncoupled equations can be specified uniquely by the initial and final conditions 
{Ro,Rt) and {ro,rt) as 

Rc{s) = RoUi{s) + RtU2{s) , 

rds) = roVi{s) +rtV2{s). (31) 

An evaluation of the path integral can be carried out which is dominated by the classical 
solutions in With these classical solutions, we write the action S[ 

S[Rc,r,] = [Mui{t)Ro + Mu2{t)Rt)rt 
- [Mu^{0)Ro + Mu2{0)Rt)ro 

+ i{au{t)rl + {au{t) + a2i{t))rort + a22{t)r^) . (32) 



1 ft ft 

(^ki{t) = - ds ds'vk{s)fL{s - s')vi{s') (33) 



where 



2 Jo Jo 

(for k,l = 1, 2) contains the effects on the system dynamics from fluctuations in the envi- 
ronment modified by the pulses. 

Using the results above, in Eq. (jl9p can be written in a compact form, 

MRt, n- 1 I i?o, ro; 0) = iV(t)e*^"— (34) 



where (a)^^ = a^^, R'^ = (i?o, Rt) and = (ro, n) 





( Uii 


U12 


)..( 




K U21 


U22 





-ni(0) ni(t) 

-U2{0) U2{t) 



(35) 
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2.3 QBM in the phase space representation 

Formally the Wigner distribution function is related to the density matrix by 

W{R, P,t) = ^ f dre-'^'pr{R + r/2, R - r/2, t). (36) 

ZTT J 

It obeys the evolution equation 

W{Ru Pu t) = J dRodPo K{Rt, Pf, t \ Ro, Pq; 0)H^(i?o, ^o, 0), (37) 
where K{R, P; t \ Rq, Pq; 0) is the propagator for the Wigner function defined by 

K{R,P;t\ Ro,Po;0) = ^J drdro e-*(^'^-^°'^°) J,(P, r; t | i?o,ro;0) . (38) 
It is given by 

K{R, P- 1 I i?o, Po; 0) = ^ / drdro e*(-^'-+^o^o+^) 

where Nw{t) = N{t) /2^J\ai\ and |a| is the determinant of a. The vector 5X = X — {X), with 

^=(f], (39) 



and 



Here S is a matrix characterizing the induced fluctuations in the system from the influence 
of the environment: 



Although many equivalent phase space representations are known |3T]. for our purpose, 
it is convenient to define the characteristic functions obtained by the Fourier transform of 
the phase space distribution functions as follows: 



XQiz,z,t) = j d'^aQ{a,a)e"'^e'"' (42) 
Xw{z,z,t)= f d'^aW{a,a)e''''e'"', (43) 



9 



where a = ^JVt/2R + i^l/2VtP . The characteristic function xq(z, for the Q-distribution 
is related to xw{.z,z), the characteristic function for the Wigner distribution, as 

XQ{z,z,t) = e ^ xw{z,z,t). (44) 

For instance, the characteristic function XQ{z-,z,t) for the harmonic QBO evolved from the 
initial ground state has the following Gaussian form: 



XQ{z,z,t) = exp 



(45) 



The time dependent coefficients are anti-normal ordered operator averages of second mo- 
ments given by (a^) = c + a, (aa^) = C + S + l/4, which have their origins in the classical 
trajectory C of a damped harmonic oscillator in Eq. ()4U|) . the induced fluctuations S from 
the environment in Eq. ()41|) . The relations of these components can be written as follows: 



8c = 


'-'22 


cf, + n'^c^,, - ^ + 2^(CiiCi2 + C21C22 


8C = 




^'c!, + ^ + CI, 


Aa = 






4S = 







(46) 



2.4 Fock states/phase space correspondence 

The harmonic QBM 

Ng 



H^^^ + H^j°°^ = na^a + v^(a + a^) ^ c„g„ (47) 

n=l 



can be viewed as a limiting case of a finite A^-level system: 

Hf^ + Hf^ = nS^S^ + (S^ + SI,)Tb (48) 

where N ^ 00. Here Tb = Yln=i^nqn and we have absorbed the factor y/2Q by defining 
Cn = \/2riCn- 5*^ is a number lowering operator whose matrix element is {S]^)ij = \/i5i^ij 
and its is a number raising operator. 

The correspondence between the Fock state representation for the pseudo spin qubits 
and the phase space representation of a density matrix is given by 

Pkiit) = [^XQ{z,z,t){k\e-'''^'e-ni). (49) 

J 71 
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^From Eq. ()49p we can directly evaluate the density matrix in the Fock representation for 
an arbitrary quantum number. The Pauli spin representation from the lowest two level Fock 
representation are obtained as 

Mt)) = Poi(t) +Pio(t) 
{(Xy{t)) = ipwit) -ipQiit) 

(a.it)) = pii(t) -poo(t). (50) 

In 13], we explained in some detail the level reduction mechanism and the relation between 
ESBM and other models of quantum computing. Various models used in the study of 
quantum computation can be obtained from our ESBM in different limits. 

At a finite temperature T, only modes in the harmonic oscillator up to ~ ksT/hQ 
are excited. At sufficiently low temperature fcsT ~ hQ, the effective excitable number of 
levels of a harmonic QBO is significantly reduced. In particular, at fc^T < HQ, we expect 
that the system is effectively reduced to two- levels. In this limit, our system and interaction 
Hamiltonians are reduced to those of the spin-boson model 

i/f + = Q{SI + ^) + S'^Tb. (51) 

Simplified model of decoherence with a pure dephasing term[25'j can be obtained as an adi- 
abatic limit of ESBM when the Brownian oscillator is frozen, i.e. f2 — with an unitary 
change of basis. This model is suitable for the study of short time dynamics but it overesti- 
mates the decoherence as we will see in the next section. Dissipative system-bath coupling 
without counter rotating terms common in many quantum optics texts 3 IJ is obtained by 
imposing a rotating wave approximation (RWA) on ESBM. 

Nb 

Hs + HsB = nStS^ + {S, + St)TB ^Hs + Hrwa = ^S^S, + dn{S,bi + Stbn),{52) 

n=l 

where bn = {uJuQu + Wn) / are bath annihilation operators. The dynamics under this 
approximation cannot capture the fast dynamics at time scales less than the natural time 
scale of the system. Thus it is inadequate in the study of generic environments. 

More generally, one can split the components of SjfS]^ as S2S2 + 5'^_2'S'7v-2) where 
Sn-2 raising and lowering operators acting on the space outside of the lowest two levels. 
Similarly, 5*^ + = S2 +S2 + S^_2 + = ^2 + ^n-2 + T2,n~2, where Sff_2 acts only on 
the space outside of the lowest two levels and {T2^N-2)nm = "\/2<5n2^m3 mixes the two sectors 
(inside and outside of the lowest two levels) and is thus responsible for the leakage. 

3 Dynamical Decoupling 
3.1 General Scheme 

Quantum error correction based on coding schemes requires many physical qubits to encode 
information [3]. In more realistic situations where only a few physical qubits are available, the 
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implementation of quantum information processing (QIP) requires more practical methods 
to prevent decoherence. This is similar to the refocusing techniques frequently used in NMR, 
where a sequence of pulses is used to modify the influence of environment on the system. 
In the bang-bang control scheme, a rapid sequence of a pair of vr pulses is used to reduce 
the effect of decoherence. In the scheme described in Sec. 2, a train of a pair of pulses with 
opposite signs to change the influence of the environment on the system. 

Pulses to diminish environmental influence In order to successfully diminish the effect 
of the environment on the system we require the pulses to satisfy the following requirements: 
First, the pulses are assumed to be strong enough so that we can ignore Hs during the pulse 
operation compared to Hp. In this strong pulse limit while keeping VLVt = {21 + 1)tt {I = 
0,1,2...) fixed, we obtain e^^p'^Hie~'^^p'^ = —Hi. Second, the pulses are assumed to be 
acting at sufficiently short intervals At. Once these conditions are satisfied, the repeated 
application of pulses will change the whole system evolution to 



where the limit with vanishing pulse intervals is taken: At ^ 0,J\f ^ oo, while t = 2AfAt 
is kept fixed. Thus under these conditions, the effect of the environment can be successfully 
eliminated. 

Note that the complete elimination of the environment is achieved only in the idealized 
limit At — > 0. In practice, it seems sufficient to add pulses faster than the environment 
characteristic scale such that the environment is effectively frozen in the duration of the 
pulses. It has been argued that this condition requires At << l/A(7y[TT], which is already 
too stringent for the known practical implementations. Thus it is worth examining whether 
this condition At -C l/Ajjv is strictly necessary. It is also important to see from our pulse 
Hamiltonian (jllj) that the total energy will be unchanged from our pulse operations. 

Pulses to minimize leakage In a multilevel system if we are only interested in limiting 
the leakage (see previous section) from the two lowest levels to the higher levels, we can 
shape the pulses to realize the decoupling operator [THj 



where is a phase factor. The condition {Ud,T2^n-2} = together with [Ud,Hs] = 



U{t = 2J\fAt) 



'^iHpT^-iAt(Hs+HE+Hi)^-iHpT^-iAtiHs+HE+Hi) 

^-iAt{Hs+HE-Hi)^-iAt{Hs+HE+Hi)___ 

^-iAt{Hs+HE-Hi)^-iAt(Hs+HE+Hi) 

^-it(Hs+HE) 



(53) 




(54) 



[Ud.He] = [Ud.S^] = Ud,S: 



yields 



U{t = 2MAt) 



Jj] ^-iAt{Hs+HE + {S^+S-f,_^+T2,N^2)TB}jJ. 



iAt{Hs+HEHS^+S'f,^2+T2,N-2)^B} 
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^ ^-iAt{Hs+HE + {S^+S%_^-T2,N-2)rB}^-i^t{Hs+HE + iSI+S^_2+T2,N-2)rB} 
^-iAt{Hs+HE + (S^+SJ,_^-T2,N^2)rB}^-iAt{Hs+HE + {S^+S^_2+T2,N-2)rB} 
_^ ^-it{Hs+HE + {SE+Sf,_2)TB} ( 

where the hmit with vanishing pulse intervals is taken: At —>■ 0, J\f —>■ oo, while t = 2MAt 
is kept fixed. The leakage can be eliminated in this limit. This method, however, will not 
change the coherence property within the computational Hilbert space. Thus decoherence 
within such space is not preventable with this method. We now discuss the dynamical 
decoupling scheme which protects the computational Hilbert space from the infiuence of the 
environment. 



3.2 Dynamical Decoupling in QBM: Quantum Control of the Un- 
certainty Relation 

We detour a bit here to show how dynamical decoupling can alter the uncertainty relation in 
quantum mechanics. The uncertainty principle ^32j is a basic requirement for any quantum 
states to satisfy. In an open quantum system, information in the system is partly lost to the 
environment, resulting in a greater uncertainty in the quantum states of the open system 
than in the original closed system. The modified uncertainty relation in QBM in Ohmic 
environment has been studied before|6j. Here we study the uncertainty relation in general 
environment. A suitable measure for the uncertainty of an open quantum system is 

= (Ai?)2(AP)2-i({A/?,AP})l (56) 

The uncertainty principle for open systems requires that for any physical state (pure or 
mixed, equilibrium or nonequilibrium) should satisfy A> 1/2. 
Using the identity 

(aat)^-(a^)(at^)=A^ + ^ + i^ + l (57) 

and (-R^) + (P^) > 1 from the usual uncertainty principle, we obtain (aa^)^ — (a^)(a'''^) > 1. 
Combining this with (ji^ and (jl^ yields the ground state population poo = ((aa"'')^ — 
(a^) (a^^))^^/^ < 1. This can also be inferred directly from a probability requirement. From 
our Fock state/phase space correspondence in (^H)), we can reconstruct the temporal phase 
space uncertainty relation arbitrarily far from equilibrium by observing a relaxation process 
in Fock space. In Fig. 1 (a), the time evolution of the uncertainty function A is plotted. We 
set M = 1 for all figures. For a state initially in a pure state, the uncertainty at t > is 
always greater than the initial value. 

Since decoupling pulses will reduce or eliminate system-environment interaction and sup- 
press the information loss in the system, the uncertainty can be used as a measure of the 
suppression of the environment- induced effect. In particular, after the application of suffi- 
ciently strong and frequent pulses, we expect that the states will obey an uncertainty relation 
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close to the initial relation, when the environment is almost completely decoupled from the 
system. 

In Fig. 1 (b), the uncertainty function at t = 1 is plotted as a function of the pulse 
frequency parameter rj = Ac/yAt/vr. Note that arbitrary pulses do not always suppress the 
uncertainty. For slow pulses, the uncertainty is increased rather than decreased. The peak 
around ~ 1 is ascribed to the resonance between the pulse and the environmental mode. For 
super-Ohmic environment, the resonance is particularly eminent, while for 1/f environment, 
the resonance is absent and the suppression of uncertainty is effective throughout the whole 
range of the pulse frequency in the figure. This phenomenon is related to so called quantum 
anti-Zeno effect. We will discuss this issue in Sec. 3.3. For pulses faster than the resonance 
r] < 1, the uncertainty function decreases as the pulses become faster and for sufficiently fast 
pulses, the initial uncertainty returns for all environments. 

3.3 Dynamical Decoupling in ESBM 

In j2ni, it was shown that the bang-bang pulses with their pulse duration close to but less 
than the threshold value tc/Auv can eliminate 1/f noises efficiently. A pure dephasing type 
model[29j was used to illustrate the case as follows. In the presence of the decoupling pulses, 
at t2j\f = 2A/'At, coherence between the two levels is given in the interaction representation 
by 

P^i(W) = e-^-(*^-^V5i(0), (58) 

where 

Dp{t2M) ~ dujl{uj) coth f ^ j tan^ (^^j (59) 

The tari^ factor contains the effect of decoupling pulses. This form indicates that decoherence 
will be suppressed for 

kuv^t < TT. (60) 

Indeed decoherence from all the modes with (4/ + l)vr/2 < cuAt < (4/ + 3)71/2 for an 
integer / will be enhanced. There is also a resonance of decoupling pulses when the pulses 
are synchronized with the particular mode in an environment causing a divergence in the 
decoherence rate that leads to rapid destruction of coherence. Bang-bang method suppresses 
1/f noise for sufficiently fast pulses, but the condition Q still needs to be satisfied. 

At a finite temperature when thermal fluctuations begin to affect the system dynamics, 
a time scale tthermai sets an additional constraint for decoupling. In order for the decoupling 
to be successful, the pulse interval needs to be smaller than the smaller value of l/Ajjv or 

tthermal • 

Here we will show that this resonance effect does not appear in more realistic ESBM 
models, t environment IS uot just a fuuctiou of Auv but dcpcuds on detailed characteristics 
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of the environment. In the case of our ESBM, t environment also sensitively depends on the 
spectral density. In more general cases, there are many more parameters to determine 
t environment, which should be determined by the spectroscopy of the environment. 

Fig. 2 shows the time evolution of the logarithm of the decay factor r{t) = —log Pi 
at T = for different environments. Sufficiently fast pulses are chosen so that the decay 
rate is suppressed by the decoupling pulses. In the super-Ohmic case, the short time decay 
is remarkable due to ultra-violet mode dominance in this environment, while in the 1// 
case, long time decay is drastic due to infrared mode abundance in this environment. The 
efficiency of decoupling directly reflects this characteristic. In the super-Ohmic case, the 
suppression of decay rate is relatively small, while in the 1// case, improvement is rather 
remarkable. 

In order to quantify how slow the pulses can be, let us introduce the parameter t] = 
AuyAt/iT. 7] < 1 means that the pulses for efficient suppression of decoherence in this 
environment need to be fast compared to the UV cutoff. This condition is difficult to be met 
experimentally. 77 > 1 means that the pulses can be reasonably slow and relatively easy to 
implement. For 1/f noise, where rj = 1.5 > 1, our DD scheme works very well. For Ohmic 
environment, 1] = 0.15 << 1, and for super-Ohmic environment, rj = 0.05 << 1, thus the 
pulses to suppress Ohmic and super-Ohmic environment need to be ultra-fast. We see that 
super-Ohmic environments can pose a serious obstacle for the implementation of coherent 
quantum operations. 

In Fig. 3 (a), coherence as a function of the pulse interval is plotted in the pure-dephasing 
case. For super-Ohmic case, the suppression is the smallest and the condition for efficient 
suppression is much more stringent than At/yAt < vr, while for 1/f noise, the pulses with At 
close to tt/Auv remain effective. However, irrespective of the environment, when the pulse 
interval gets longer, A^/yAt ~ vr, there is a crossover from decoherence suppression to deco- 
herence accentuation. There is a sharp resonance point at A^yAt = vr, where decoherence 
rate diverges for all environment and the condition A^/yAt < n need be strictly enforced. 

In Fig. 3 (b), the effect of decoupling pulses in Eq. (fTUj) is studied in our ESBM. The 
main difference from Fig. 2 (a) is that: (i) there is no sharp resonance as seen in Fig. 3 (a) in 
any environment although the crossover from decoherence suppression to accentuation can 
still be seen, (ii) the condition A^yAt < vr is not necessary for the suppression of 1// noise. 
Even very slow pulses AjyyAt >> vr can suppress it. The second condition (ii) expressed 
before in Eq. (2) is particularly important in that it is often too stringent for practical 
implementation with qubits. Note that the crossover behavior expressed in (i) is already 
manifest in the uncertainty relation plotted in Fig. 1(b). Similar effect also appears in the 
study of quantum Zeno effect where sufficiently frequent measurements suppress the popu- 
lation decay[33 . When the measurement is performed less frequently, there is a crossover 
from quantum Zeno effect to the so called quantum anti-Zeno effect, where measurement 
accelerates the decay^Q. 
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3.4 Strong, long and slow pulse limit 

We now study the limit of long (r oo) and slow (At oo) pulses, opposite to the cases 
considered so far. In this limit, the pulse Hamiltonian can be written as the interaction 
between the system Hamiltonian and the constant field: 

Hp HsV. (61) 

Let us consider the time evolution operator acting on the eigenstate of the total Hamiltonian: 

U{t)\N) =e-"'^'\N), (62) 

where H\N) = En\N). We take the large V limit and consider the first order correction 
with respect to V^^. Using the nondegenerate perturbation theory, we see that |A^) = 
\n) + y-i|^)(™r-rerfion)^ ^ y^^ ^ ^^|^^ + Hj + HeIu) + 0{V-^), whcrc |n) and e„ = fin 
are eigenstates and eigenvalues of the Hamiltonian Hs- By expanding ()62|) around these 
unperturbed states, we obtain in the large V limit, 

U{t)\N) = e-'^'^'ln) + 0{V-') 

(63) 

for all A^. This implies 

U{t) = e-^«^+i)^«+^^)* + O(V^-i) 

(64) 

Thus in the limit of large V, the system-bath interaction Hj is eliminated. This effect resem- 
bles the quantum Zeno effect [HS] in which the decay is suppressed by frequent measurements. 
Indeed a parallel argument with the quantum Zeno effect case is possible (I64p and one can 
show that there is no relaxation in this limit [21]. From (jMj) 

Pnn{t) = {n\p{t)\n) 



Tr 



Uit)mU-\t)\n){n 



= Tr[/)(0)|n)(n|] = P„(0). (65) 

The difference here from the quantum Zeno effect in [33j caused by repeated measurements 
is that the time evolution of the system in this case is unitary. The effect of the external field 
in (jnH) is to slow down the decay process. It can in principle provide an alternative method 
to dynamical decoupling although V is often uncontrollable as it is fixed by the design 
and fluctuations of V induce an additional source of decoherence. In Fig. 4, the excited 
state population is plotted as a function of the strength of the constant external field V. 
For Ohmic and super-Ohmic cases, frequency renormalization (Ohmic) and frequency and 
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mass renormalization (super-Ohmic) terms are taken into account by introducing counter 
terms [H^ before applying V. These terms dominant at ultra-short times are not suppressed 
by the constant field in (j6T| while they can be suppressed by the pulse decoupling (fTO|) as we 
saw in Sec. 3. 2 and 3.3. In Fig. 4, the decay can be suppressed from coupling to constant large 
external field V after the renormalization although overall suppression is not so remarkable 
as decoupling by pulses . In pulse decoupling, the effect of the environmental mode with 
frequency u « 1/At is suppressed, while in the constant field case, for any finite V, all the 
environmental modes upto u; ~ 1/t contribute to decoherence. 

3.5 Conclusion 

In this paper, we propose a dynamical decoupling scheme for error prevention or deterrence 
based on the effective spin-boson model we introduced recently. A train of decoupling pulses 
with alternating signs creates a fictitious instantaneous Hamiltonian evolution. Our analysis 
does not make any of the conventional approximations such as the Born, Born-Markov, RWA 
or adiabatic approximation, thus enabling us to probe the system dynamics in a much wider 
parameter regime including the low temperature, ultra-short time scale and strong external 
field conditions which are likely to be more relevant for quantum information processing. We 
studied this method in a general environment with Ohmic, super-Ohmic and 1// spectral 
density functions. In all cases, in our decoupling scheme, there is no saturation of decoherence 
rate due to resonance as reported in previous studies using bang-bang controls. We saw the 
smooth crossover from suppression to accentuation of decoherence in our model that extends 
in a wide range of pulse durations. In particular, our method can suppress 1// noise with 
much slower pulses than previously studied. Due to the shielding of high energy modes in the 
environment as a result of defect motion, decoherence caused by telegraph noise is weaker 
at short times than that by 1// noise from the harmonic environment discussed in this 
paper. Thus we expect our pulse decoupling method could also provide efficient suppression 
of telegraph type noises. 

For efficient suppression of decoherence from an Ohmic environment , the pulse parameter 
needs to satisfy rj < 1 and the pulses need to be fast compared to the UV cutoff, which is 
difficult to implement. For 1/f noise, r] > 1, thus the pulses can be reasonably slow and our 
scheme is relatively easy to implement. However, for super-Ohmic environment, ?7 << 1 and 
the pulses need to be ultra-fast. Thus the most serious obstacle for the implementation of 
quantum computing may be the presence of super-Ohmic environments. 

We also studied the effect of strong constant field on the population decay. The decay 
process appears to slow down but the suppression is not so evident compared to the fast 
pulse decoupling. 

In reality, we may not know enough information about the environment to identify the 
possible sources of decoherence [22] • But one can cleverly reverse the argument here and 
propose to use our decoupling method to obtain information about the environment spectrum 
by measuring the decoherence rate with changing pulse intervals. For instance, if the value of 
7] at which the crossover between decoherence suppression and accentuation occurs is greater 
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than 1, the environment has sub-Ohmic spectrum. If smaller than 1, super-Ohmic. 

At this point, we are not aware of any experimental results on super-Ohmic environments. 
Since an electromagnetic field seen by a charged electron due to charge-field interaction in 
a three dimensional space is a super-Ohmic environment (HH] we expect this type of environ- 
ment to be rather ubiquitous than peculiar. While the 1// noise due to charge-fluctuations 
in Coulomb forces are frequently observed and more dominant, super-Ohmic noise due to 
charge-photon interaction is likely to coexist with 1// noise. However, the short time dom- 
inant characteristics of super-Ohmic noise will make the actual detection very difficult. 

In quantum computation with multiqubits, the entanglement among qubits is used as a 
resource for the speed-up of computation based on quantum algorithm. In the presence of 
entangling operation, the propagation of error in the whole computing space is unavoidable. 
To overcome this, fault tolerance is required in quantum error correction schemes jSZ]- The 
spreading of non-Markovian noise, in particular, noise of the super-Ohmic type in the com- 
putation architecture can be a serious threat to many proposed fault-tolerant computation 
schemes [38 . 

Our results are readily applicable to the study of quantum aspects in nanomechanical 
resonators [5^ for which the system is commonly treated as a harmonic oscillator. ESBM[3] 
captures the essential short time dynamics common in many qubit models, in which their 
computational Hilbert space is obtained by the lowest two levels of the underlying multi-level 
structure. It also allows us to calculate the precise non-Markovian time evolution of leakage. 
The leakage is commonly neglected or crudely estimated in most of the popular theoretical 
models. Another application of our work is to superconducting Josephson junction qubits 
I1T1I12], in particular, the phase qubit mo dels [IB]. 
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Figure 1: (a) The exact temporal evolution of the uncertainty function ()5(i|l in units of 
of a system initially in its ground state interacting with an Ohmic and super-Ohmic 
environment is plotted. The dotted line indicates the minimum uncertainty at 1/4. All 
the curves start off from this value at t = 0[ns] and stay above the line. The effect of 
vacuum fluctuations of the environment sets in at t ~ 0.1 [ns] and the effect of thermal 
fluctuations of the environment at t ~ 1/T = 0.3 — 0.5[ns] for Ohmic case. Both effects 
set in earlier for super-Ohmic case. Other parameters are u = l[GHz], 7 = 0.1[GHz], 
Auv = lO[GHz]. (b) Same system now in the presence of decoupling pulses at t = 0.25[r;,s]. 
Other parameters are n = l[GHz]. T = U.5[mK], 7 = 0.1[GHz], Auv = lOO[GHz] for 
Ohmic, T = U.5[mK], 7 = 0m02[GHz], Auv = 50[GHz] for super-Ohmic, T = 0[mK], 
7 = 0.1[GHz], Auv = lOO[GHz] for 1/f environments. 
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Figure 2: Plots of the logarithm of the decay factor at T = 0, Q = l[GHz]. Solid (Dashed) 
curves are decay with (without) decoupling pulses. Parameter values are: Ajjv = lOO[GHz], 
7 = 0.l[GHz], N = 100 for Ohmic, and Auv = lOO[GHz], 7 = 0.5[GHz], N ^ 30 for 1/f, 
Auv = 3Q[GHz], 7 = QM[GHz], N = 100 for super-Ohmic case. 



23 




Figure 3: (a) Coherence as a function of tlie pulse interval at T = is plotted, t = 0.25[ns]. 
Parameters are: Auv = 100[GHz], Air = 1. 7 = 0.2[GHz] for 1/f, 7 = 0.125[GHz] 
for Ohmic, and 7 = 0.005[GHz] for super-Ohmic case, (b) Plots of the first excited state 
population at T = and t = 0.15[ns]. Other parameter values are: Auv = 100[GHz], 
n = 10[GHz], 7 = 0.1[GHz] for Ohmic, and Auv = lOO[GHz], n = 15[GHz], 7 = 0.5[GHz] 
for 1/f, Auv = 50[GHz], fl = 15[GHz], 7 = 0.01[GHz] for super-Ohmic case. Parameter 
values are chosen to fit all the curves in the parameter range in the figure. For the same 
parameter values, while their qualitative behavior remains similar, the differences in their 
values are more striking. 
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Figure 4: Excited state population as a function of V is plotted at t = 0.15[ns]. The initial 
state is the first excited state. Auv — 100[GHz]. Q = l[GHz\. Parameter values are: 
7 = 0.1[GHz] for Ohmic, and 7 = 0.1[GHz] for 1/f, 7 = 0.001[GHz] for super-Ohmic case. 
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